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{ ! ABSTRACT 

Context. There are numerous extrasolar giant planets which orbit close to their central stars. These 'hot-Jupiters' probably formed in the outer, 
I cooler regions of their protoplanetary disks, and migrated inward to ~ 0. 1 AU. Since these giant planets must have migrated through their inner 

■ systems at an early time, it is uncertain whether they could have formed or retained terrestrial planets. 

T^j- Aims. We present a series of calculations aimed at examining how an inner system of planetesimals/protoplanets, undergoing terrestrial planet 
7-H , formation, evolves under the influence of a giant planet undergoing inward type II migration through the region bounded between 5-0.1 AU. 
CO ■ Methods. We have previously simulated the effect of gas giant planet migration on an inner system protoplanet/planetesimal disk using a N- 
body code which included gas drag and a prescribed migration rate. We update our calculations here with an improved model that incorporates 
a viscously evolving gas disk, annular gap and inner-cavity formation due to the gravitational field of the giant planet, and self-consistent 
evolution of the giant's orbit. 

Results. We find that > 60% of the solids disk survives by being scattered by the giant planet into external orbits. Planetesimals are scattered 
pIZh ' outward almost as efficiently as protoplanets, resulting in the regeneration of a solids disk where dynamical friction is strong and terrestrial 
Oh' planet formation is able to resume. A simulation that was extended for a few Myr after the migration of the giant planet halted at 0.1 AU, 
q resulted in an apparently stable planet of ~ 2 m e forming in the habitable zone. Migration-induced mixing of volatile-rich material from 

■ beyond the 'snowline' into the inner disk regions means that terrestrial planets that form there are likely to be water-rich. 

"T^ Conclusions. We predict that hot- Jupiter systems are likely to harbor water-abundant terrestrial planets in their habitable zones. These planets 
C$ • may be detected by future planet search missions. 
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1. Introduction. Initially, it was thought reasonable that the inward migra- 
tion of a giant planet would be so disruptive of the material 
it passed through as to clear the swept zone completely, pre- 

About one quarter of the extrasolar planetary systems discov- cludin 8 the formation of any inner system planets. However, 

ered to date contain a so-called 'hot-Jupiter' - a gas giant planet thls Vlew was a conservative assumption, often made in the 

orbiting within 0.1 AU of the central star (e.g. Butler et al. su PP ort of speculative astrobiological arguments (e.g. Ward & 

2006). It is improbable that these planets formed within such Br°wnlee 2000; Lineweaver 2001; Lineweaver et al. 2004). 

a hot region of the protoplanetary disk, and is most hkely that Agreement as to the outcome failed to materialize from the 

they originated further out beyond the nebula snowline and first modeling studies of the process, the conclusions of which 

moved inward via type II migration (e.g Lin & Papaloizou vaned from the occurrence of terrestrial planets in hot-Jupiter 

1986; Ward 1997; Bryden et al. 1999; Nelson et al. 2000). s y stems bein 8 m 8 nl y unllkel y ( Arm itage 2003), through pos- 

Since both giant planet formation and type II migration re- sible but rare (Mandell & Sigurdsson 2003), to commonplace 

quire the nebular gas to still be present, these processes are con- (Raymond et al. 2005a). None of these studies, however, actu- 

strained to occur within the first few million years of the disk aU y simulated terrestrial planet formation simultaneously with 

lifetime (Haisch et al. 2001), well within the ~ 30 - 100 Myr giant planet migration. Their disagreement about the likely out- 

timescale thought to be required to complete the growth of come can be traced t0 assumptions made about the timing of 

terrestrial planets (e.g. Chambers 2001; Kleine et al. 2002; giant planet formation and migration. 

Halliday 2004). Thus, these migrating giant planets must have The first study to model inner system planetary accre- 

traversed the inner system, including its habitable zone, before tion in the presence of a migrating giant was that of Fogg & 

any planet formation there was complete, raising the question Nelson (2005) (hereafter referred to as Paper I). This work used 

of what effect such a disturbance would have on the growth of N-body simulations to examine oligarchic and giant-impact 

terrestrial-like planets. growth (Kokubo & Ida 2000) in a protoplanet/planetesimal disk 
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based on the minimum mass solar nebula model of Hayashi 
(1981), extending between 0.4 - 4 AU. Five scenarios were 
considered, corresponding to five different ages for the inner 
planet forming disk at the point when a giant planet was as- 
sumed to form at 5 AU and migrate in to 0. 1 AU. 

In all five of their scenarios, Fogg & Nelson (2005) found 
that the majority of the disk solids survived the passage of the 
giant planet, either by being shepherded inward of the giant, 
or by being scattered by the giant into excited exterior orbits. 
This partition of solid material was shown to vary with the level 
of dissipative forces present, which decline with disk matu- 
rity, favoring shepherding at early times and scattering at late 
times. Within the portion of the disk compacted inside the in- 
creasingly restricted volume interior to the giant, accretion was 
found to speed up greatly at late times resulting typically in a 
~ 3 — 15 M ffi planet forming inside 0.1 AU. The similarity of 
these objects to the recently identified class of 'hot-Neptune' 
planets (McArthur et al. 2004; Butler et al. 2004; Santos et al. 
2004b; Vogt et al. 2005; Rivera et al. 2005; Bonfils et al. 2005; 
Udry et al. 2006) was noted and discussed. The fate of the ma- 
terial scattered into external orbits was not subjected to fur- 
ther calculation, but it was noted that, although the remaining 
solids surface density was reduced from pre-migration-episode 
values, ample material remained to provide for the eventual ac- 
cretion of a set of external terrestrial planets, including within 
each system's habitable zone. Fogg & Nelson (2005) therefore 
concluded that the assumption that hot-Jupiter systems are de- 
void of inner system terrestrial planets is probably incorrect, 
and that planet formation and retention both interior and exte- 
rior to the hot-Jupiter is possible. We note that similar results, 
relating to the formation of planets interior to a migrating giant, 
have been reported by Zhou et al. (2005). 

In this paper, we extend the model introduced in Paper I by 
improving the realism with which gas dynamics is simulated. 
A 1-D evolving viscous gas disk model is linked to the N-body 
code that: 1) allows the gas to deplete over time via viscous ac- 
cretion onto the central star; 2) allows an annular gap to form 
in the vicinity of the giant planet; 3) includes the creation of 
an inner cavity due to dissipation of propagating spiral waves 
excited by the giant planet; 4) self-consistently drives the giant 
inward. Compared to the unevolving gas disk assumed in Paper 
I, this new model reduces the strength of dissipation present 
in all scenarios, especially in regions close to the central star 
and the giant. We examine and re-interpret the fate of the disk 
solids under these changed circumstances and look at whether 
our hypothesis of hot-Neptune formation remains robust. We 
also examine the post-migration evolution of the outer scat- 
tered disk of solids, and find that terrestrial planets do form in 
the habitable zone. Another issue we examine is the extent to 
which volatile-rich matter, originating from beyond the snow- 
line, is driven into the inner system and is mixed into the sur- 
viving planet-forming material. We find that substantial mixing 
occurs, such that any terrestrial planets that form are likely to 
be water-rich bodies hosting deep, global oceans. 

The plan of the paper is as follows. In Section 2 we outline 
the additions to our model and the initial conditions of the sim- 
ulations; in Section 3 the results are presented and discussed; 



in Section 4 we consider some caveats and future model im- 
provements, and in Section 5 we offer our conclusions. 



2. Description of the model. 

We model planetary accretion using the Mercury 6 hybrid- 
symplectic integrator (Chambers 1999), run as an N + N' sim- 
ulation, where we have N protoplanets embedded in a disk of 
N' "super-planetesimals" - particles that represent an idealized 
ensemble of a much larger number of real planetesimals. The 
protoplanets (and the giant when it is introduced) interact grav- 
itationally and can accrete with all the other bodies in the sim- 
ulation, whereas the super-planetesimal population is non-self- 
interacting. These latter objects however are subject to a drag 
force from their motion relative to the nebular gas. A detailed 
outline of these aspects of our model is given in Sect. 2.1 & 
2.2 of Paper I and we continue here to discuss the additional 
features we have introduced. 

2.1. Improved nebula model. 

In Paper I we assumed a steady state gas disk model with a 
constant surface density profile oc r~ L5 . The migration rate of 
the giant was prescribed from a calculation of the local viscous 
evolution timescale. More realistically, the quantity of nebular 
gas should decline due to viscous evolution and accretion onto 
the central star, progressively depleting the inner disk. The gas 
density should also decrease in the vicinity of the giant due to 
the generation and dissipation of density waves at Lindblad res- 
onance positions, clearing an annular gap in a zone where the 
planetary tidal torques dominate the intrinsic viscous torques 
of the disk. A consistent calculation of the type II migration 
rate should involve the back-reaction to these tidal torques. 

To account for these processes we model the gas disk 
by solving numerically the disk viscous diffusion equation 
(Pringle 1981), with modifications included to account for the 
tidal influence of an embedded giant planet. Such a method has 
been used previously in studies that attempt to explain the sta- 
tistical distribution of exoplanetary orbits through type II mi- 
gration and disk dispersal (Trilling et al. 1998; Armitage et al. 
2002; Alibert et al. 2005). The simplest technique for includ- 
ing the effect of the planet is the impulse approximation of Lin 
& Papaloizou (1986), where wave dissipation is assumed to 
occur close to the planet. A more sophisticated treatment of 
the problem is the WKB approximation (Takeuchi et al. 1996) 
which involves summing the torque contributions from a se- 
ries of Lindblad resonances in the disk. The former technique 
was adopted in the studies cited above as it requires consid- 
erably less computation and generates comparable results. We 
follow this approach here, but in order to include the effect of 
non-local dissipation of waves that travel far into the disk we 
also include the WKB approximation torques due to the waves 
launched at the two innermost and outermost Lindblad reso- 
nances. 
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We assume a MMSN-type protoplanetary disk around a 
M, = 1 M Q star (Hayashi 1981). The initial surface density 
of solids is: 



£ s (r) =/neb/ice£l (y^jj) 



(1) 



where r is radial distance from the central star, / ne t is a nebular 
mass scaling factor, ~L\ = 7 g crrT 2 and the ice condensation 
coefficient /; ce = 1 for a < 2.7 AU (the distance chosen for the 
nebula 'snowline') and /j ce = 4.2 for a > 2.7 AU. As in Paper 
I, we set / neb = 3. 

The initial surface density of gas is: 



Eg(r) = /neb/gas^l (j^jj) 



(2) 



where / gas is the gas to dust ratio which we take to be / gas = 
240. 

Given a nebular radial temperature profile of r neb = 
280(r/ALTr 1/2 , the sound speed of a solar composition gas in 
cgs units is: 



c s = 9.9 x 10 



/ ^neb Y 

\280K/ ' 
and the gas scale height is: 



h = 0.047 



(l Au) 



AU . 



(3) 



(4) 



Since the kinematic viscosity in an alpha-disk model is v = 
ah 2 £l, where Q. is the local Keplerian angular velocity, we take 
the turbulent shear viscosity function of the disk to be: 



v = 9.84x 10 

and the viscous evolution time: 

2 f^ 2 



r v =-(£) (aQ)-^lO-VV, 48^^) 



(5) 



(6) 



In all the models presented here we assume a disk alpha vis- 
cosity of a = 2 x 10 -3 , giving t v (5 AU) ~ 120 000 years. 
We solve the diffusion equation for S g (r) in the form: 



5Zg 
dt 



ld_ 

r dr 



id, i \ 2AI„ri 
3ri — (v2 g r5j L_ 



dr 



(GM»)5 3n(GM„)-- 



(7) 



where the first term in square brackets describes the diffusion of 
gas under the action of internal viscous torques (Pringle 198 1); 
the second term describes the impulse approximation of the lo- 
cal tidal interaction of the planet with the disk, with A being the 
specific torque exerted by the planet (Lin & Papaloizou 1986); 
and the third term, which derives from the WKB approximation 
(Takeuchi et al. 1996), is included to account for more distant 
angular momentum transfer via the damping of waves launched 
from the innermost and outermost two Lindblad resonances, 
with T being a summation of the torque densities exerted by 
these waves. Because their launch sites stand off a substantial 



distance from the planet, these waves are expected to be lin- 
ear and not damped locally in the disk. Their angular momen- 
tum content is therefore deposited in the disk through viscous 
damping as they propagate. 

The exchange of angular momentum between the planet 
and disk leads to a radial migration of the planet at a rate: 



da _ / a \ J 1 
dt ~ \GM,j m p 



An 



rAI g dr + 2 



Tdr 



(8) 



where m p is the mass of the planet, a is its semi-major axis 
and rin and r out are the inner and outer boundaries of the disk 
respectively. 

The rate of specific angular momentum transfer to the disk 
in the impulse approximation is given by Lin & Papaloizou 
(1986) as: 



2 GM,( r 



(9) 



where q = m p /M„, and A p = max(h, \r - a\). 

The total torque density exerted on the disk in the WKB 
approximation via the damping of waves excited by the planet 
is (Takeuchi et al. 1996): 



T(r) = J] T m (r) , 



(10) 



where m is the mode number of the mth order Lindblad reso- 
nance at: 



1 



a , 



(11) 



where use of the minus sign gives the radial distances of 
the inner resonances r IL and the plus sign those of the outer 
resonances r 0L . Since only the innermost and outermost two 
Lindblad resonances are accounted for here, we take m - 2,3 
for the resonance positions interior and m = 1, 2 for those ex- 
terior to the planet. 

The torque density is calculated from the radial gradient of 
the angular momentum flux F m (r): 



TJr) 



dF m (r) 
dr 



(12) 



which is given in Takeuchi et al. (1996) as: 



F m (r) = F mQ exp 



"I! 



m(Q p - Q) 



p) 2 ) 



kdT 



(13) 



where F m (r) = for rrL < r < roh, C is the bulk viscosity (set 
here to zero), k = Q for a Keplerian disk, Q p is the angular 
velocity of the planet and k(r) is the radial wavenumber: 



k{r) 



m 2 {Q. - Q p ) 2 - k 2 



(14) 
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Takeuchi et al. (1996) give this approximation for the angular 
momentum flux originating at a given mth order resonance: 



Table 1. Data describing initial solids disk set-up 



F na = - m %^r L ){— 



I Gm 



2K (^) + K 1 ^j_^„ u (2±l) 



(15) 



where Ko and Ki are modified Bessel functions, t) m j is the 
Kronecker delta function, and the upper component of the + 
is used for the inner and the lower component for the outer res- 
onances respectively. 

The parameter f c is a torque cutoff function (Artymowicz 
1993) given by: 



fc = 



1 



H(l +4£ 2 ) 



2//K (2///3) + Ki(2///3) 



2K (2/3) + K!(2/3) 



(16) 



where H = (1 + f) 1/2 and £ = m(c s /aQ p ). 

The evolution of the nebular gas is computed by solving 
Eq. 7 with an explicit finite-difference technique on a grid with 
a cell width Ar oc yjr. The resulting type II migration forces 
on the giant planet are computed from Eq. 8 by deriving an 
instantaneous time scale, a I a, and inserting this into Eq. 6 of 
Paper I. Strong eccentricity damping for the giant planet is as- 
sumed, with the damping time scale being l/50th of the radial 
migration timescale. As in Paper I, we have neglected the ef- 
fects of type I migration and associated eccentricity damping 
(Ward 1997; Tanaka et al. 2002; Tanaka & Ward 2004). 

The gas disk adopted here extends from an inner radius 
ri n = 0.0025 AU to an outer radius of r out = 33 AU with 
an initial 2 g (r) profile given by Eq. 2. Since we are consid- 
ering a 3 x MMSN disk, / ne b = 3 and the initial disk mass is 
M gas = 0.0398 M « 42 Mj. Note that this initial disk mass is 
greater than the 2 x MMSN of gas assumed in Paper I; however 
since that amount was kept fixed, the gas present in the simula- 
tions presented here falls below this level after ~ 140 000 years. 
The boundary conditions for the computation were £ g (r in ) = 0, 
representing gas accretion onto the central star and at r out the 
radial velocity of the gas was set to zero. 

To correctly couple the evolving gas disk algorithm with the 
N-body code, synchronization of their respective time-steps is 
necessary. In each simulation sub-run (see Sect. 2.3) the sym- 
plectic N-body time-step T n body was fixed whereas the gas disk 
time-step T gas is adaptive and taken to be: 



T g as = nun 



w Ar(i) 
M0I 



(17) 



where ;' is the grid cell label and w is a coefficient of order unity 
that is tuned to ensure computational stability. Including all the 
torques given above, the gas radial velocity is: 



A + 



Sgr 1 / 2 dr 



vV 1/2 ] 



(18) 



Thus, if T gas > * 1 



nbody then we set r gas = T nbo d y ; if T g 



^"nbody 

then the gas disk is evolved for INT(T n body/i" gas ) steps of dura- 
tion Tg as plus an extra step of MOD(T n body/i" gas )- 





Rocky Zone 


Icy Zone 


Total 




0.4-2.7 AU 


2.7-4.0 AU 


0.4-4.0 AU 


M solld 


9.99 M ffi 


14.8 M e 


24.8 M e 


^proto 


0.025 M e 


0.1 Me 




N 


66 


9 


75 




0.0025 Me 


0.01 Me 




N' 


3336 


1392 


4278 


/proto 


0.17 


0.06 


0.1 



2.2. Radial mixing of solid material. 

At early times, the solid component of a young protoplane- 
tary nebula will exhibit a radial pattern of chemical compo- 
sition, controlled by the temperature-dependent condensation 
sequence of a variety of metals, rock minerals and ices. As the 
planetary system grows and evolves, phenomena such as dy- 
namical spreading, gas drag induced orbital decay and resonant 
interactions can cause a radial mixing of material. According to 
one school of thought (e.g. Morbidelli et al. 2000) the original 
matter that condensed in the Earth's orbit is thought to have 
been dry (for an alternative opinion see Drake & Righter 2002) 
and the origin of the Earth's water, and its D/H ratio, can be 
explained if ~ 10% of the planet's mass is composed of car- 
bonaceous chondrite-type material, originating from between 
2.5 - 4 AU, and ~ 10% of the water gained thereby is retained 
at the end of accretion. N-body simulations of terrestrial planet 
formation from disks that extend out close to the orbit of Jupiter 
are supportive of this idea and all demonstrate substantial mix- 
ing of water rich material into the inner disk (e.g. Morbidelli et 
al. 2000; Chambers 2001; Raymond et al. 2004, 2005a,b). 

In Paper I it was noted that one consequence of the inward 
migration of a giant planet is the shepherding of planetesimals 
that are damped by gas drag (Tanaka & Ida 1999), and the trap- 
ping of bodies at first order mean motion resonances. Hence the 
outer, more volatile-rich, regions of the protoplanetary disk are 
actively mixed into its inner regions. However, the composition 
of planetesimals and protoplanets, and their accretion products, 
were not logged in our simulations. 

We amend this here by labeling all particles with a compo- 
sition based on its original position in the disk, and summing 
the composition of protoplanets as they grow. In a similar man- 
ner to most other studies, we assume a crude three-phase initial 
radial composition with rocky material originating at < 2 AU, 
material characteristic of chondritic meteorites between 2 - 
2.7 AU, and trans-snowline material at > 2.7 AU. For conve- 
nience, these phases are referred to as "dry", "damp" and "wet" 
respectively but, in contrast to other studies, we do not assign 
an actual water mass fraction to them. Instead we use our re- 
sults to make more generalized observations and predictions. 

2.3. Initial conditions and running of the simulations. 

In common with Paper I, we assume a nominal age for the 
protoplanetary disk of 0.5 Myr, this being the t — start 
date for our simulations, and a mass of three times that of 
the MMSN. Our reasons for choosing a more massive nebula 
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Table 2. Overall solids disk data: after 0.1-1.5 Myr of evolution 




250000 500000 750000 1000000 1250000 1500000 

t(years) 



Fig. 1. Evolution of the mass of the gas disk. The mass of gas (in 
Jupiter masses) remaining at the launch point for each of the five mi- 
gration scenarios is indicated. 




r(AU) 



Fig. 2. Evolution of the gas surface density within the inner 10 AU of 
our simulated disk. The upper solid line is the r~ l 5 X g -profile for a 
3 x MMSN disk. The lower solid curves, in descending order, are the 
profiles at 0. 1 , 0.25, 0.5, 1 .0 and 1 .5 Myr respectively. The dashed line 
is the fixed Z g -profile assumed in Paper I. 



stem from the fact that core-accretion theories of giant planet 
formation require an enhanced density of solid material in or- 
der to grow a critical core mass before the nebular gas is lost 
(Lissauer 1987; Thommes et al. 2003) which might be sup- 
ported by the observation that hot-Jupiters are found predom- 
inantly around stars more metal-rich than the Sun (Santos et 
al. 2004a; Fischer & Valenti 2005). As runs proceed, the initial 
gas disk described by Eq. 2 evolves according to the algorithm 
outlined in Sect. 2.1, initially without the presence of a giant 
planet (i.e. A = 0, T = in Eq. 7). 

Simultaneously, Mercury 6 evolves an initial disk of solid 
material, extending radially from 0.4 - 4.0 AU. We use the 
same two-component initial solids disk as was used in Paper 
I, generated according to the profile of Eq. 1 and in line with 



Time (Myr) 


0.1 


0.25 


0.5 


1.0 


1.5 


Scenario ID 


I 


II 


III 


IV 


V 


M soM (M e ) 


24.8 


24.8 


24.8 


24.8 


24.8 


m max (M e ) 


0.29 


0.37 


0.63 


1.40 


2.13 


N 


50 


47 


44 


38 


34 


N' 


4031 


3665 


3036 


2342 


1964 


/proto 


0.19 


0.26 


0.37 


0.51 


0.59 



the oligarchic growth picture of Kokubo & Ida (2000), which 
we assume to be a reasonable description of the state of the in- 
ner disk at 0.5 Myr. Data for this initial disk model are shown 
in Table 1 which gives, for zones interior and exterior to the 
snowline, values for the total mass of solid material M so i;d, the 
number and mass of protoplanets N and m proto , and the num- 
ber and mass of super-planetesimals N' and m s _ p i. Note that as 
super-planetesimals act as statistical tracers for a much more 
numerous population of real planetesimals, their mass is much 
greater than that of a real planetesimal, but for the purposes of 
calculating gas drag, each super-planetesimal is treated as if it 
is a single 10 km radius object of realistic mass. The parameter 
/proto. at the foot of Table 1, is the mass fraction of the solids 
disk contained in protoplanets and we use this here as a rough 
measure of the evolution of the disk, taking / pro to = 0.5 to de- 
note the transition between oligarchic and giant impact growth 
regimes. 

From t — 0, we run our new model for 0.1, 0.25, 0.5, 1.0 
and 1 .5 Myr in the absence of the giant, with T n b d y = 8 days 
and a simulation inner edge of r* =0.1 AU. (Note that r, does 
not necessarily denote the physical radius of the central star, but 
it is the distance interior to which particles are removed from 
the simulation and their masses added to the central star.) The 
evolution of the nebular mass over this time span and the par- 
ticular nebular mass at each of these five epochs are shown in 
Fig. 1. The gas surface density profiles resulting at these times 
are shown in Fig. 2. It is apparent that gas drains onto the cen- 
tral star very rapidly at first, as the density gradient relaxes to a 
shallower profile. Compared to the unevolved profile, order of 
magnitude reductions in gas density occur within the disk's in- 
ner regions. Data for the evolved solid components are given in 
Table 2 and include the values of m max , the mass of the largest 
protoplanet to have evolved in each case. (In contrast to Paper I 
we do not run our model to 3.0 Myr, as there remains too little 
gas at this time to provide for our giant planet's envelope.) The 
advance of planetary growth with time is indicated by the pro- 
gressive increase of m max and / proto and the reduction in particle 
numbers as super-planetesimals are accreted and protoplanets 
merge. When compared to the equivalent data from Paper I, we 
note that no planetesimals are lost beyond the inner edge of the 
simulation because of much lower gas densities near the central 
star. 

The five type II migration scenarios studied here are con- 
structed from the five evolved disks indicated in Fig. 1 and 
summarized in Table 2. Scenarios I - III take place whilst the 
solids disk remains in its oligarchic growth phase (/ pro to < 0.5) 
whereas Scenarios IV - V have just entered the final giant im- 
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Scenario I 

Scenario II 

Scenario III 

Scenario IV 

Scenario V 



t (years) 

Fig. 3. Semi-major axis evolution of the giant planet in each scenario 
from the launch time (the top row in Table 2) to the time at which it 
arrives at 0. 1 AU 

pact stage of growth (f ?mto > 0.5). A giant planet of mass 
0.5 M } is placed into each simulation at 5.0 AU after removing 
0.4 Mj of gas from between 3-7 AU. The giant then proceeds 
to clear out an annular gap in the disk and migrates inward 
according to the method described in Sect. 2.1. The runs are 
halted once the giant reaches 0. 1 AU (see Fig. 3). For Scenarios 
I - HI, this takes ~ 1 10 000 - 120 000 years, close to the pre- 
diction of Eq. 6. The process takes longer to complete in the 
cases of Scenarios IV - V (~ 150000 and 220000 years re- 
spectively) because by the time of the appearance of the giant 
planet the gas disk is substantially depleted and is less effective 
at driving migration. In order to better model processes when 
the giant migrates down to small radial distances, we contract 
the simulation inner edge down to a realistic T-Tauri star radius: 
r, = 0.014 AU = 3 R Q . The initial timestep chosen for the sym- 
plectic integrator was T nbody = 8 days, but it was necessary to 
reduce this at late times as material is driven into closer orbits. 
Hence each scenario was divided into a number of sequential 
sub-runs with T n b d y being adjusted at each re-start so as to keep 
the timestep close to one tenth the orbital period of the inner- 
most object. Since planetesimals in this new model suffer less 
orbital decay due to gas drag, it was possible to conduct these 
runs with a higher value of r n body than in Paper I. However, this 
advantage was negated at late times as the adaptive r gas falls 
steeply as the giant planet moves within 1 AU. From t — 0, the 
scenarios presented here each required 4-8 weeks of 2.8 GHz 
CPU-time for completion. 

3. Results. 

We begin describing the results by focussing on the evolution 
during the migration of the giant planet, discussing one specific 
case in detail before examining differences between the various 
runs. We then go onto describe the evolution of Scenario V 
after giant planet migration has halted, focussing on the issue of 



terrestrial planet formation in the scattered disk, and the likely 
composition of planets that form there. 

3.1. Evolution during giant planet migration 

The results of all scenarios showed a number of behavioral fea- 
tures in common. As in Paper I, we illustrate these first by de- 
scribing the results of Scenario I in detail. We then proceed to 
examine how the results differ between scenarios (dependence 
on disk maturity) and how the results differ from those of Paper 
I (dependence on an evolving gas disk). 

3.1 .1 . Typical features of a run. 

The typical effects of a migrating giant planet on an inner solids 
disk observed from Paper I were as follows: 1) shepherding 
of planetesimals; 2) capture of objects at first order mean mo- 
tion resonances; 3) acceleration of accretion interior to the gi- 
ant with possible hot-Neptune formation; and 4) generation of 
a scattered exterior disk. To a greater or lesser extent, these 
features were also observed in our new simulations. Four snap- 
shots of the evolution of Scenario I are illustrated in Figs. 4-7 
showing the mass, inclination and eccentricity of objects, and 
the gas surface density vs. semi-major axis. The original prove- 
nance of the protoplanets (interior or exterior to the snowline) is 
denoted by the shading of its symbol as described in the caption 
to Fig. 4. In the case of a merger between rocky and icy proto- 
planets, this shading is determined by that of the most massive 
of the pair. 

An early stage in the evolution of Scenario I, 20 000 years 
after the introduction of the giant planet, is shown in Fig. 4. The 
giant has opened a ~ 0.75 AU gap in the gas and has migrated 
inward to 4.37 AU, shepherding the outer disk edge at the 4:3 
resonance, now at 3.61 AU. Capture of objects at the 3:2 and 
2:1 resonances, at 3.33 and 2.75 AU respectively, is apparent 
from eccentricity spikes visible in the upper panel and a cluster- 
ing of protoplanets in the lower panel. Even at this early phase, 
before the giant has entered the original confines of the inte- 
rior disk, one protoplanet and a handful of super-planetesimals 
have been scattered into external orbits. 

The system midway through the run, 60 000 years after the 
introduction of the giant planet, is shown in Fig. 5. The giant 
has now migrated to 2.72 AU, putting the positions of the 3:2 
and 2:1 resonances at 2.07 and 1.71 AU respectively. Strong ex- 
citation of protoplanetary orbits is now apparent at these loca- 
tions, as is the build-up of scattered material in external orbits. 
The primary mechanism of this expulsion is evident from the 
behavior of material captured at resonances. Continuous reso- 
nant pumping results in orbits becoming eccentric enough to 
eventually intersect the orbit of the giant planet. A series of 
close encounters with the giant then follows, eventually result- 
ing in a final encounter where the object is accreted or expelled 
into a non-intersecting exterior orbit. 

An advanced stage of Scenario I, 100000 years after the 
introduction of the giant planet is shown in Fig. 6. The giant 
planet is now at 0.70 AU and the 3:2 and 2:1 resonances are at 
0.54 and 0.44 AU respectively. A substantial scattered external 
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Fig. 4. Scenario I at 20 000 years after the start of giant planet migration, showing the mass, inclination and eccentricity of objects. Small black 
dots represent super-planetesimals; white filled circles are rocky protoplanets; grey filled circles are icy protoplanets and the large highlighted 
grey filled circle is the giant. Objects plotted between the dotted lines in the upper panel have orbits that intersect the orbit of the giant. The 
location of the 2:1, 3:2 and 4:3 resonances with the giant are indicated. Gas surface density is read on the right hand axis of the lower panel, 
the upper grey curve being the unevolved profile at t = and the lower black curve being the current profile. 



disk has now formed and sufficient gas remains in this early 
scenario to rapidly damp the orbits of scattered planetesimals. 
An impression of the scattering process in action is given by 
the five protoplanets currently crossing the giant's orbit. The 
interior disk is compacted to high surface densities, but now 
that strong first order resonances with the giant are influential 
throughout its width, and gas densities have fallen by a factor 
of > 10, there is a noticeable dynamical stirring of its entire 
remaining population. Nevertheless, accretion has speeded up 
in this shepherded zone with the growth of one protoplanet of 
1.17 me at 0.53 AU outrunning that of its neighbors. 

Scenario I is terminated at 114000 years after the start of 
migration when the giant planet arrives at 0.1 AU and the sys- 
tem is illustrated at this point in Fig. 7. Two thirds of the orig- 
inal solids disk mass has survived the migration episode - the 
great majority of this residing in the scattered exterior disk. 
Most of the remainder has been accreted by the giant planet. 
Just ~ 4% of the original disk mass remains interior to the gi- 
ant and none of this is in the form of large bodies. The rapidly 
accreting interior protoplanet seen in Fig. 6 continued its in- 



ward progress close to the 3:2 resonance position with its or- 
bit being well-damped by strong collisional damping and dy- 
namical friction from planetesimals and smaller protoplanets 
(see Fig. 10 in Sect. 3.1.4). However, ~ 2000 years before its 
demise, the protoplanet drifted outward and became captured 
at the 4:3 resonance. At this location, both accretion and dy- 
namical friction were reduced allowing the protoplanet's orbit 
to become progressively more eccentric. At 109220 years the 
protoplanet, now weighing in at 2.41 m e , collided with and 
was accreted by the giant. The five other less massive interior 
planets visible in Fig. 6 grew very little, remaining between 
~ 0.03 - 0.2m e , and in due course one of them impacted the 
giant and four were scattered. 

To emphasize the above description of Scenario I, the sur- 
face density evolution of the disk and its accretion rate are 
shown in Fig. 8. The left hand panel shows the disk surface den- 
sity profile (obtained by summing all protoplanets and super- 
planetesimals in 0.1 AU width bins) at 20000, 40000, 60000, 
80000 and 100000 years after the start of migration; the right 
hand panel plots the amount of mass accreted onto protoplan- 
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Fig. 5. Scenario I at 60 000 years after the start of giant planet migration. The giant has now moved inward to 2.72 AU. Increasing excitation of 
the orbits of protoplanets captured at resonances is apparent, as is the build-up of matter scattered into external orbits. 



ets only (including protoplanet mergers) every 1 000 years for 
the duration of the run. In the 2 S plot, two surface density en- 
hancements are clearly visible as spikes at the 3:2 and 2:1 res- 
onances and are seen to grow whilst moving inward. At 80 000 
years, these have almost merged into one: the shepherded por- 
tion of the original disk having by now been squeezed into a 
dense ring. By 100000 years, most of this mass is now con- 
fined within 0.6 AU and E s here has risen to ~ 500 g ctrT 2 
which is off the vertical scale in the figure. This amounts to an 
increase by a factor of ~ 10 over the previous, undisturbed, disk 
surface density, but is only about half the increase seen in the 
equivalent Scenario presented in Paper I. The effect of this disk 
compaction process is visible in the accretion rate plot. Mass 
accretion rises significantly after 80 000 years due to both the 
high values of E s and the fact that much of this mass now re- 
sides in a zone where dynamical times are shorter. However, 
the large, terminal, accretion rate spike described in Paper I is 
not reproduced here (compare Fig. 8 with Fig. 6 of Paper I). 
This is because close to the end of that previous simulation 
a 15.65 m e hot-Neptune was assembled in a dramatic phase of 
runaway accretion interior to 0. 1 AU. In the case presented here 
the compacted interior disk is only half as dense and is much 



less well damped (note the large difference between the upper 
and lower gas density curves in Fig. 7) and, whilst a proto- 
planet does grow to 2.41 m ffl in this region, as described above, 
it does not survive and is accreted by the giant planet. We de- 
scribe the formation and fate of interior planets in more detail 
in Sect. 3.1.4. 

3.1 .2. Dependence on the maturity of the inner disk. 

The reason for running five scenarios through a progressively 
more mature inner disk is to see if the timing of migration has 
any systematic effect on the results. This is possible as the par- 
titioning of the solids disk between inner and outer remnants is 
influenced by the level of damping that particles are subject to, 
which declines with age. In Paper I, where E g (r) is fixed, this 
occurs as a side effect of accretion: as planetesimals are ac- 
creted by protoplanets (f ploto increases), fewer small particles 
remain that are subject to gas drag and which can exert dynam- 
ical friction. In these latest simulations, since we now have an 
evolving gas disk, the strength of gas drag on susceptible parti- 
cles also declines with time and is particularly marked close to 
the central star and giant planet. 
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Fig. 6. Scenario I at 100000 years after the start of giant planet migration. The giant has now moved inward to 0.70 AU. Five protoplanets are 
currently crossing the orbit of the giant. The scattered disk has grown and a > 1 me planet is accreting within the compacted interior disk. 

Table 3. Fate of the disk mass at the end of each run. 



Scenario 


I 


II 


III 


IV 


V 


Total Initial Solids (M e ) 


24.81 


24.81 


24.81 


24.81 


24.81 


Total Surviving Solids (Me) 


16.60 (67%) 


16.69 (67%) 


17.40 (70%) 


21.23 (86%) 


20.22 (81%) 


Interior Surviving Solids (M ffi ) 


0.88 (4%) 


0.65 (2%) 


1.00 (4%) 


0.84 (3%) 


0.31 (1%) 


N , ,/proto 


0,0 


0,0 


0,0 


0,0 


0,0 


Exterior Surviving Solids (M©) 


15.72 (63%) 


16.04 (65%) 


16.40 (66%) 


20.39 (82%) 


19.90 (80%) 


fpioto 


39, 0.27 


29, 0.28 


33, 0.42 


31,0.63 


23, 0.66 


Accreted by Star (M e ) 


0.01 (0.04%) 


0.01 (0.04%) 


0.11 (0.4%) 


0.0 (0%) 


0.0 (0%) 


Accreted by Giant (M ffi ) 


8.20 (33%) 


7.85 (32%) 


6.77 (27%) 


3.41 (14%) 


4.59 (19%) 


Ejected (M e ) 


0.00 (0%) 


0.26(1%) 


0.51 (2%) 


0.17(1%) 


0.0 (0%) 



Data describing the fate of the solids disk mass at the end 
of each scenario are shown in Table 3. Disk mass that is lost is 
either accreted by the central star, ejected from the system, or 
accreted by the giant planet; that which survives is partitioned 
between bodies orbiting interior or exterior to the final orbit of 
the giant planet at 0. 1 AU. 

In all scenarios, a negligible quantity of mass was ejected 
or lost to the central star. However, a significant fraction of 
the disk mass (14 - 33%) was accreted by the giant, espe- 
cially towards the end of the migration. At these late times, 



planetesimals are shepherded into the partially evacuated in- 
ner regions of the gas disk where gas drag is less effective at 
damping orbital perturbations from the giant and growing pro- 
toplanets. Once a planetesimal strays into the annular gap in 
the gas containing the giant, gas drag vanishes and accretion 
or scattering by the giant follows. This increased excitation of 
the shepherded planetesimal population, and the thinning down 
of their number, renders them less effective at damping proto- 
planetary orbits via dynamical friction and collisions. Hence, at 
late times the orbits of the remaining interior protoplanets also 
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Fig. 7. Scenario I at 114 000 years after the start of giant planet migration. The giant planet has migrated to 0. 1 AU. Most interior mass has 
been lost after the most massive interior protoplanet impacts the giant. 63% of the original solids disk mass now resides in exterior orbits. 



tend to destabilize, with one of the same two fates in store. A 
trend can be seen in Table 3 for the giant planet to accrete less 
material with disk maturity. This occurs because as the disk 
ages the gas density and the solids mass fraction in small bod- 
ies both decline, resulting in less dynamical dissipation of both 
planetesimals and protoplanets. Less matter is shepherded in 
such mature disks so there is less of an interior remnant for the 
giant to accrete from at late times. 

In all scenarios, a large majority of the disk solids are found 
to survive the migration episode - over two thirds of the orig- 
inal inventory. Table 3 shows that there is essentially no trend 
with disk maturity in the partitioning of surviving mass be- 
tween interior and exterior remnants. Just a few percent of the 
mass remains interior to the giant in all cases. When the giant 
planet migrates through a disk undergoing oligarchic growth 
(/proto < 0.5; Scenarios I - III) ~ 65% of the original disk mass 
survives by being scattered into the exterior disk. This frac- 
tion increases for disks undergoing giant impact-style growth 
(/proto > 0.5; Scenarios IV - V) to > 80%, not because more 
mass remains at < 0. 1 AU but because less mass is accreted by 
the giant at late times. 



3.1 .3. Dependence on an evolving gas disk. 

The salient dynamical behaviors of solids disk particles such 
as shepherding, resonant capture, scattering by the giant planet 
and eventual partition into interior and exterior remnants are 
observed generally in the results of both the present model and 
those of Paper I. However, the introduction of an evolving gas 
disk causes the relative predominance of these outcomes to dif- 
fer. This is because both the principal sources of dissipation, 
dynamical friction and gas drag, fall with time, whereas in the 
model from Paper I only the former declines. There are five 
systematic differences between the results of Paper I and those 
presented here. 



1 . Much less mass is lost to the central star. 

2. Much more mass is accreted by the giant planet. 

3 . Disk partitioning between inner and outer remnants is much 
less sensitive to disk maturity. 

4. The protoplanet mass fraction (/ pro to) m the exterior scat- 
tered disk is lower. 

5. Interior hot-Neptune-type planets grow to smaller masses 
and do not survive. 
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Fig. 8. Surface density evolution (left hand panel) and accretion rates (right hand panel) for Scenario I. Growing surface density peaks at the 
2:1 and 3:2 resonances sweep through the inner system ahead of the giant. Accretion rates increase after ~ 80000 years within the compacted 
portion of the disk. 



The reasons for the first two items are the large reductions 
in gas density close to the central star and in the vicinity of the 
giant planet. Planetesimals do not have time to spiral into the 
central star and feel little eccentricity damping when close to 
the giant because of the low gas drag in these regions. The rea- 
sons for the last three items stem from the fact that the system is 
less dissipative so scattering behavior predominates over shep- 
herding at all epochs. Interior remnants are consistently much 
smaller than outer remnants which always include > 60% of 
the original solids disk mass. The increased tendency for plan- 
etesimals to scatter means that there is less of a fractionation 
of planetesimals from protoplanets between interior and exte- 
rior remnants. More planetesimals are expelled into the exterior 
disk and hence its overall mass fraction contained in protoplan- 
ets (/ptoto) is lower. 

These tendencies are illustrated in Fig. 9 where the percent 
of the original solids disk mass surviving at the end of the sim- 
ulations presented here are compared with the simulation re- 
sults from Paper I. Interior and exterior remnants are plotted as 
dashed and solid lines respectively; gray lines are the results of 



Paper I and black lines are those of the present work. It can be 
seen that partitioning varies strongly with disk maturity in the 
case of a steady-state gas model and weakly, if at all, when the 
nebula gas is allowed to viscously evolve. Since gas drag does 
not decline in the former case, the shepherding of planetesimals 
is more influential and more mass remains in the interior rem- 
nant in early scenarios, which is mostly contained in surviving 
hot-Neptune-type planets. Late scenarios from Paper I behave 
more similarly to the ones presented here as a greater fraction 
of the solids disk mass is contained within protoplanets, which 
do not feel the gas, and their source of dissipation, dynamical 
friction, is weaker due to a decline in super-planetesimal num- 
bers. 1 

The increased tendency in the present model for planetes- 
imals to scatter into the exterior disk can also be appreciated 
from Fig. 9 by comparing, for the two models, the ratio of 
the protoplanet mass fraction in the final exterior disk to that 



1 Note that Scenario V is less comparable to its counterpart in Paper 
I as it is initiated at a different time, 1.5 Myr earlier. 
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Fig. 9. Comparison of the results of Paper I (grey lines and x sym- 
bols) and the results of this paper (black lines and + symbols). Dashed 
lines and solid lines are the percentage of original disk mass found in 
the interior and exterior remnants respectively. The ratios of the pro- 
toplanet mass fraction in the final exterior disk to that of the original 
disk, / pr() , (final) // proto (initial), are indicated by the x and + symbols 
and are read off the right hand y-axis. 



of the original disk, / pro , (nnal)//p roto (initial). If the scatter- 
ing process does not discriminate between protoplanets and 
planetesimals, and there is no further accretion after scatter- 
ing, this ratio would remain at / proto (final)//p roto (initial) = 1. 
Subsequent accretion will raise this ratio, but not by much 
over the timescale considered, especially as accretion rates are 
reduced in the exterior disk due to its dynamically excited 
state. Reading from the right hand y-axis, Fig. 9 shows that 
/proto (final)/ /p roto (initial) < 1.5 for all of the present model re- 
sults whereas much higher values are obtained from Paper I. 
Protoplanets are preferentially scattered into the exterior disk 
when there is strong and unvarying gas drag, whereas fraction- 
ation is much less marked when gas drag declines over time 
and falls to near zero in the vicinity of the giant planet. 

3.1.4. Formation and survival of hot-Neptunes 

A striking feature of the results of Paper I was the growth and 
survival, interior to the final orbit of the giant, of hot-Neptune 
and hot-Earth type planets ranging between ~ 2 - 16 m®. No 
such planets are found to survive in the runs of our present 
model. Accelerated accretion in the shepherded disk is ob- 
served with some protoplanets growing to ~ 1-3% (see 
Figs. 6 and 8), but in each scenario their eventual fate was to 
impact the giant planet just before the end of the simulation. 
After formation, these planets become locked in a mean motion 
resonance (3:2 or 4:3) with the giant, and subsequent eccen- 
tricity pumping during migration leads to collision and merger 
with the giant shortly before migration halts. However, had we 
halted giant planet migration earlier, at ~ 0.25 - 0.5 AU in- 
stead of 0.1 AU, these simulations would predict the existence 
of hot-Earths (e.g. see Fig. 10). 



Fig. 10. Blow-up of the inner 1 AU of Scenario I, showing eccentricity 
vs semi-major axis 105 000 years after the start of giant planet migra- 
tion. Protoplanetary masses are indicated in m ffi and the locations of 
the 2:1 and 3:2 resonances are shown with arrows. The giant planet is 
at 0.48 AU and our hot-Earth candidate is a 1.93 m e planet at 0.34 AU 
with e ss 0.03. Five lower mass protoplanets are in the process of being 
scattered into the exterior disk. (Compare with Fig. 6 which shows the 
situation 5 000 years earlier.) As described in Sect. 3. 1 . 1 , the hot-Earth 
candidate is eventually accreted by the giant planet ~ 4 000 years later 
once it has migrated further inward to ~ 0.26 AU. 



We have performed additional calculations, very similar to 
those presented here, in which we switched off the effect of dis- 
sipating spiral waves travelling into the disk (i.e. T(r) = in 
Eq. 7). This term helps create an inner cavity in the gas disk 
as the planet migrates inward. These calculations produced in- 
terior planets ranging in mass ~ 2 - 4 ra e during giant planet 
migration. A variety of final outcomes were noted: 1) some of 
them were accreted by the giant planet; 2) others were scat- 
tered externally into stable orbits at a ~ 0.4 AU; and 3) a few 
remained in interior orbits, typically close to 0.076 AU (the 3:2 
resonance). These results, taken alongside the five scenarios al- 
ready presented in this paper and those from Paper I show a dis- 
tinct trend: a strongly dissipating gas disk interior to the giant 
planet leads to the formation of fairly massive, hot-Neptune- 
like planets which survive; a gas disk of lower density leads to 
the formation of lower mass interior planets that often do not 
survive. 

We note at this point that our models currently neglect some 
potentially important sources of dissipation due to general un- 
certainty about how planetary formation can proceed in their 
presence, such as type I migration (Ward 1997; Tanaka et al. 
2002; Tanaka & Ward 2004) and the circularization of orbits 
due to stellar tides. Associated with the former is strong eccen- 
tricity and inclination damping of low mass planets which may 
facilitate the survival of inner planets. Simulations are currently 
underway to examine this possibility. Tidal damping of orbits 
however is unlikely to have a significant effect on our results 
because all shepherded protoplanets are scattered or accreted 
by the giant at a > 0.1 AU. Whilst tidal circularization times 
at these radii are uncertain, due to the uncertainty of Q, they 
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probably range from ~ 10 8 - 10 10 yr (Goldreich & Soter 1966), 
orders of magnitude greater than the millenia it takes for our 
giant planets to traverse the final 0.5 AU of their migration. 

3.1 .5. The exterior scattered disk. 

In Paper I it was noted that the exterior solids disks generated 
by the giant planet migration were dynamically excited, de- 
pleted of planetesimals, and spread over a greater radial extent 
than the original disk. Whilst subsequent planet formation still 
seemed possible in such a disk, it was predicted to occur over 
significantly extended timescales. However, in the Paper I sce- 
narios only ~ 1 - 2m ffi of planetesimals were scattered into the 
exterior disk, whereas in the present model this quantity rises 
to ~ 8 - Urn®. In addition, planetesimals are not scattered as 
widely as protoplanets and their excited orbits damp rapidly 
when remote from the influence of the giant planet. Whilst the 
decline in gas drag with time lessens this trend in later scenar- 
ios, in all cases the inner regions of the scattered disk remain 
well populated. 

There is a greater similarity between the current models in 
the outcome of scattering of the protoplanetary population and 
data for the external protoplanets are shown in Table 4, giving 
their number, mean and maximum masses and orbital inclina- 
tions, and their mean, minimum and maximum semi major axes 
and eccentricities. As expected, protoplanetary numbers fall 
and masses rise with disk maturity, an effect largely due to prior 
accretion before the appearance of the giant planet. Allowing 
for stochastic events, such as giant impacts and strong scat- 
terings, these data are much the same in Paper I, as are the 
mean semi-major axes and eccentricities (a « 5 AU, e « 0.5). 
However the minimum values of the semi-major axes and ec- 
centricities (a m i n and e m ; n ) are lower than in Paper I. This is be- 
cause late-shepherded protoplanets tend to scatter rather than 
assembling into hot-Neptunes and their resultant orbits damp 
more rapidly as many more planetesimals are available to exert 
dynamical friction. 

3.2. Post-migration terrestrial planet formation 

What do the results of the five scenarios presented here have 
to say about the probablity of forming terrestrial planets in the 
scattered disks ? Considering the Scenarios I - V overall, if the 
factor of long-term importance is the mass distribution, rather 
than initial dynamics, then such planets should form and await 
discovery in hot-Jupiter systems. This point is made in Fig. 1 1 
where the total solids mass with semi-major axes between 0.75 
- 1.75 AU, before and after giant planet migration, is plotted 
for each scenario. A clear trend is visible for less matter to be 
found in this region with increasing disk maturity as it is more 
widely scattered. Mass dispersal in late scenarios, however, is 
partially offset by the pre-existing inward evolution of mate- 
rial in more mature disks which enhances the mass present in 
inner regions. It is also offset by the fact that fairly massive pro- 
toplanets are scattered into the external disk in late Scenarios 
(IV and V - see table 4), so that terrestrial planet formation 
in the scattered disk has already received a significant boost 




Scenario 



Fig. 11. The total solids mass between 0.75 - 1.75 AU, before and 
after the giant planet migration, plotted for each scenario. 



from accretion prior to and during migration. In all cases, more 
than 2 m e of planet forming material remains in the 'maximum 
greenhouse' habitable zone (Kasting et al. 1993) of the system 
after the passage of the giant planet. 

Extending our runs for the additional 50 - 100 Myr it 
would take to form a completed external planetary system is 
beyond our current computational capabilities. The principal 
difficulty is the presence of the hot-Jupiter at 0.1 AU which 
limits Tnbody to an excessively low value. We have, however, ex- 
tended Scenario V for 2 Myr and find that (apparently stable) 
terrestrial planets do form in the habitable zone of this system. 

As the gas density is low when migration of the giant planet 
halts in Scenario V, we assume rapid removal of the remnant 
gas via photoevaporation and evolve the system in the absence 
of gas. Fig. 12 shows the region between 0-4 AU for the 
scattered disk of Scenario V at the point when the giant planet 
reaches 0.1 AU (top panel). Protoplanets of substantial size 
have been scattered into the external disk. The result after a fur- 
ther 2 Myr of gas-free accretion is shown in the bottom panel. 
Dynamical friction is sufficient to cause general damping of 
protoplanetary orbits at the expense of excitation of the plan- 
etesimals. A significant accretion event has been a giant im- 
pact and merger between the 1 .42 & 0.73 m e protoplanets at 
~1.5 AU resulting in a 2.15 % body at a = 1.47 AU with 
e = 0.13. This planet, which lies within the 'maximum green- 
house' habitable zone, no longer crosses the orbit of any of its 
neighbors and is probably a long-term survivor. There has been 
some rearrangement of the intersecting protoplanetary orbits 
beyond 2 AU with two protoplanets moving outward and cir- 
cularizing at a > 4 AU, off the right hand side of the figure, 
to compensate for the 0.79 AU inward movement of a massive 
2.2 m e planet, now at a = 2.25 AU with e = 0.14. This planet 
is currently crossing the orbit of a smaller 1.07 m ffi body, the 
long-term fate of which will probably be accretion by one of 
its two larger neighbors or incremental scattering by them to 
a safe distance. The emergence of a stable terrestrial planetary 
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Table 4. Data describing the external surviving protoplanets at the end of giant planet migration. 



Scenario 


N 


"Vote (Me) 


Wmax (M e ) 


a(AU) 


"mm (AU) 


"max (AU) 


e 




^max 


i ° 


'max 


I 


39 


0.11 


0.25 


5.43 


0.45 


37.92 


0.51 


0.15 


0.99 


8.85 


40.45 


II 


29 


0.16 


0.74 


4.99 


0.71 


15.25 


0.49 


0.045 


0.91 


8.59 
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Fig. 12. Eccentricities of bodies within 4 AU at the end of Scenario V (top panel) and after a further 2 Myr of gas-free accretion (bottom panel). 
Protoplanets are shown as white circles and are labeled with their mass in trig,. Super-planetesimals are indicated by black dots and the grey 
blob at 0.1 AU is the hot- Jupiter in its post-migration orbit. Gas density at the end of Scenario V is shown as the black curve in the top panel 
(read off the right-hand axis), the grey curve being the unevolved profile at t = 



system from such a configuration seems highly probable, and 
this result provides a clear prediction that terrestrial planets will 
be found in the habitable zones of hot-Jupiter systems. 

What will be the physical nature of such planets ? The is- 
sue of disk and protoplanet composition after migration is dis- 
cussed in detail in Sect. 3.3. Here we note, however, that the 
2.15 and 2.20 m ffl planets observed in Fig. 12 are composed 
of ~ 30% and 45% of material originating from beyond the 
snowline, respectively. Assuming that trans-snowline planetes- 
imals and protoplanets contain 75% water, a naive prediction is 
that these planets will contain between 20 - 30 % water by 
mass. The current water inventory of the Earth is estimated to 



be about 0.05 % by mass. Evidently terrestrial planets form- 
ing in the habitable zones of hot-Jupiter systems are likely 
to host deep global oceans - essentially being 'water-worlds' 
(Kuchner 2003; Leger et al. 2004; Raymond et al. 2006), even 
if significant loss of volatiles occurs during high-impact accre- 
tion. 

This prediction about water content, however, depends on 
the giant planet forming out beyond the inner edge of the snow- 
line, rather than at its inner edge. The closer to the inner edge 
the giant forms, the smaller the amount of volatile-rich mate- 
rial that will be shepherded inward, and the lower the degree of 
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"damp" 

"wet" 




Fig. 13. Composition of the original solids disk in Scenario I (top 
panel) compared with the composition of the scattered disk generated 
through giant planet migration (bottom panel). The key is explained 
in the text. 



little less smooth in Scenarios IV & V which have had more 
time to accrete substantial protoplanets from local material be- 
fore the appearance of the giant planet. 

In contrast to other studies of water delivery to terrestrial 
planets (e.g. Morbidelli et al. 2000; Chambers 2001; Raymond 
et al. 2004, 2005a,b, 2006), we do not assign an actual water 
mass fraction to our compositional phases. This is because our 
simulations do not extend to the completion of planetary accre- 
tion and our MMSN-type original solids surface density profile, 
which includes a large step-increase at 2.7 AU (see Eq. 1), is 
different from those adopted in the above-cited papers making 
a detailed comparison difficult. Typically, we find that at the 
end point for early-migration scenarios (I - III) originally dry 
protoplanets that have found their way into the scattered disk 
contain - 25% of trans-snowline material, but are surrounded 
by a large quantity of volatile-rich planetesimals from which 
to accrete further. In later scenarios (IV - V), where the sur- 
viving protoplanets are more mature, this range increases to 5 
- 70%, and whilst there is less remaining mass in small bod- 
ies to sweep up in these systems, the accretion of volatile-rich 
material is not yet complete. 

In general we can predict that the final terrestrial planets 
that form in hot-Jupiter systems are likely to be much more 
volatile-rich than the Earth, those in the habitable zone ending 
up with > 10% of their material originating from beyond the 
snowline. Accretion in the external disk, however, may pass 
through an early high velocity phase before completion with 
some protoplanets potentially losing much of their volatile in- 
ventory during giant impacts (Asphaug et al. 2006; Canup & 
Pierazzo 2006). Since water-rich planetesimals are abundant 
throughout our external disks, however, protoplanets that have 
lost their original volatile endowment in a catastrophic event 
should be able to re-acquire some water before planetary ac- 
cretion is complete. Dry terrestrial planets appear improbable 
in hot-Jupiter systems whereas Earth-like worlds and planets 
with deep global oceans (Kuchner 2003; Leger et al. 2004; 
Raymond et al. 2006) may well be commonplace. 



volatile enrichment experienced by terrestrial planets forming 
after migration of the giant. 

3.3. Migration-driven compositional mixing. 

A result found in all Scenarios I-V is homogenization of the 
solids disk composition. This mixing occurs as the giant planet 
drives material inward from the outer part of the disk whilst 
generating the external disk via random scattering. An example 
is given in Fig. 13 which shows, for Scenario I, the composition 
of the solids disk before migration (top panel) and that of the 
surviving external disk after migration (bottom panel). In each 
case the total mass of solid material is plotted in 0.5 AU width 
bins with the histograms labeled "dry", "damp" and "wet" rep- 
resenting rocky material originating at < 2 AU, material char- 
acteristic of chondritic meteorites between 2 - 2.7 AU, and 
trans-snowline material at > 2.7 AU, respectively. It is evi- 
dent that a large amount of material from beyond the snowline 
is shepherded into the inner regions before being left behind. 
Compositional mixing is similar in other scenarios, although a 



4. Discussion. 

The models presented here and in Paper I have only explored 
a small region of the parameter space relevant to this problem, 
and have inevitably adopted assumptions that simplify or omit 
potentially important physical processes. We discuss some of 
these issues below and their possible implications for our re- 
sults. 

(/) Giant planet mass & radius. In all simulations we have 
adopted a giant planet mass of 0.5 Mj, with a density of 
1 g ctrT 3 and a radius of 0.85 Rj. Whilst the giant planet can 
accrete solid matter as it migrates inward, it does not accumu- 
late any more gas and so its mass changes by < 5% during a 
run. Hot-Jupiters, however, come in a variety of masses and al- 
though their average m sin i is less than a Jupiter mass, more 
massive examples are known. We have not yet tested the effect 
of increasing the mass of our migrating giant planet, but expect 
that this would result in an enhanced and more widespread scat- 
tering of material into the external disk, and more accretion by 
the giant. 
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We have assumed a radius for our giant planet that is represen- 
tative of a fully contracted state, which may not be realistic so 
soon after its formation. Giant planets of approximately Jovian 
mass contract to radii ~ 2 - 3 Rj during their rapid gas ac- 
cretion phase (Papaloizou & Nelson 2005), before cooling and 
contracting toward the Jovian radius over longer time scales. 
Adoption of a larger radius would probably result in greater 
accretion of solids by the giant, with gas-drag-sensitive plan- 
etesimals being preferentially accreted. 

(ii) Migration halting mechanisms. The reason why hot- 
Jupiters cease their inward migration close to the central star 
and come to rest at < 0. 1 AU is unknown. Various mechanisms 
have been suggested, such as the giant planet moving into a 
central magnetospheric cavity in the gas disk (Lin et al. 1996), 
or fortuitous dispersal of the gas disk due to photoevaporation 
(Trilling et al. 1998; Armitage et al. 2002; Alibert et al. 2005). 
We do not rely on any particular mechanism here and merely 
note that if the former operates then any of our scenarios may 
be appropriate, whereas if the latter is true then our latest sce- 
narios (IV & V) are more likely to be realistic as they operate 
at a time when most of the gas disk has already been lost. In 
either case, we expect an external disk to be generated during 
giant planet migration and renewed terrestrial planet formation 
to follow. The dynamical state of this scattered disk, however, 
will be affected by photoevaporation of the gas disk as the gi- 
ant migrates. The effect of this on terrestrial planet formation 
during and after giant planet migration will be examined in a 
future paper. 

(Hi) Type I migration. We have not included the effect of type 
I migration (Ward 1997; Tanaka et al. 2002; Tanaka & Ward 
2004), which would operate most effectively on sub-gap open- 
ing bodies of > lm ffi , in any of our published models to date and 
have assumed that it does not play a major role. However, we 
discussed in detail its possible influence in Paper I and confine 
ourselves here to the further comment that even if type I migra- 
tion operates with a reduced efficiency to that expected from 
the canonical model, which seems likely, it might still have a 
significant effect because of the sensitivity of our model's out- 
come to the dissipation exerted on solid bodies. Type I interac- 
tions would exert damping forces on protoplanets additional to 
dynamical friction. This could have the effect of reducing the 
eccentricity excitation of protoplanets captured at resonances 
with the giant planet, hence reducing their tendency to intersect 
the orbit of the giant and be scattered into the external disk. The 
overall effect might be to increase the shepherding effect of the 
giant planet at the expense of scattering, enhancing the frac- 
tion of surviving material found interior to the giant planet at 
the end of the simulation. If this material avoids migration into 
the central star, it is possible that inclusion of type I migration 
forces will assist in the formation and survival of hot-Neptune 
and/or hot-Earth systems. We are currently running simulations 
to examine this possibility. 

(iv) Planetesimal size evolution. For computational simplic- 
ity we have assumed a uniform planetesimal population with 
radii of 10 km. In reality there would be a distribution of plan- 
etesimal sizes determined from a balance between their rates 
of accumulation and fragmentation. In regions of the disk that 
are dynamically cold, the mean planetesimal size would grow 



via binary mergers, whilst in regions that are more dynami- 
cally excited destructive collisions could result in the planetes- 
imal population being ground down into smaller fragments. 
Fragmentation of the planetesimal population in the context 
of the oligarchic growth regime however is not necessarily 
an obstacle to planet formation and may actually assist plan- 
etary growth by supplying protoplanets with a more strongly 
damped feedstock, enhancing the effect of gravitational fo- 
cussing (Chambers 2006). A particularly challenging environ- 
ment for planetesimal survival in our simulations in in the com- 
pacted region of the disk between the 4:3 and 2:1 resonances 
with the giant planet. In this region, planetesimal surface den- 
sities are enhanced (see the left hand panel in Fig. 8) and the 
population as a whole is strongly stirred generating eccentric- 
ities as high as 0.2 - 0.3 (see Figs. 5 and 6). Random veloc- 
ities of several km per second, far in excess of planetesimal 
escape velocities, are indicated and mutual planetesimal colli- 
sions would result in fragmentation. What this means for our 
model is unclear, but it is possible that the fragment popula- 
tion, which would be much more strongly affected by gas drag, 
would evolve inward rapidly, thereby escaping the dynamically 
excited zone. It could then be gathered efficiently by protoplan- 
ets and accreted. We speculate that the overall effect of this pro- 
cess might be to reduce the planetesimal fraction of the external 
disk material, but not to reduce its overall mass as protoplanets 
should be scattered there with the same, or even greater, effi- 
ciency. 

(v) Eccentric giant planet. In this paper we have assumed 
that the giant planet is on a high circular orbit as it migrates. 
Nagasawa et al. (2005) have shown that the interior disk can be 
pre-stirred through the action of a sweeping secular resonance 
if the giant is on a modestly eccentric orbit. This would modify 
the accretion history of the inner disk. 

5. Conclusions. 

In this paper we have presented the results of simulations that 
model terrestrial planet formation during and after the migra- 
tion of a gas-giant planet to form a 'hot-Jupiter'. This work is 
an extension of our previous work (Fogg & Nelson 2005), with 
improvements being made by modelling the viscous evolution 
of the gas disk, gap-formation and inner cavity formation in 
the gas disk due to the gravitational influence of the planet, and 
self-consistent type II migration. A popular belief has been that 
hot-Jupiter systems are unlikely to host terrestrial planets, as 
migration of the giant planet through the terrestrial planet zone 
was expected to sweep that region of planet forming material. 
We find, however, that the majority of this mass survives the 
migration episode as an interior or exterior disk remnant from 
which terrestrial planet formation can resume. This occurs via 
a combination of shepherding of the original solids disk ahead 
of the giant, and random scattering of the majority of this com- 
pacted material into orbits external to the giant. The net effect is 
not a disappearance of planetary building blocks from the inner 
system, but rather a stirring and mixing of material originally 
formed at different radial distances. 

Now that more realistic gas dynamics have been included in 
our model, generating partial cavity formation close to the cen- 
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tral star, gap formation in the vicinity of the giant planet's orbit, 
consistent type II migration rates, and a decline in the overall 
mass of gas present with time, we have found the following 
qualitative differences between the results presented here and 
in our previous work: 

1 . Scattering is favored over shepherding irrespective of the 
maturity of the inner solids disk at the epoch of giant planet 
migration. This occurs because of the rapid loss of gas from 
the disk interior to the giant, and in the vicinity of the giant 
planet's orbit, which reduces the damping of planetesimal 
trajectories. The efficiency of shepherding and dynamical 
friction is therefore reduced, increasing the probability that 
a given body will come close enough to the giant to be scat- 
tered into an external orbit. In all our scenarios > 60% of 
the original solids disk material survives in a regenerated 
external disk. 

2. The principal mass loss mechanism is accretion by the gi- 
ant planet, rather than the central star. This is because gas 
densities close to the central star are reduced by 2 - 3 orders 
of magnitude, suppressing the gas-drag-induced orbital de- 
cay of planetesimals and allowing the giant planet to catch 
up with and sweep through even the innermost solids disk 
material. Most accretion by the giant occurs at late times 
when scattering by the giant becomes less effective due to 
the system being contained deep within the star's gravita- 
tional potential well. 

3. Hot-Neptune and/or hot-Earth survival is not favored in our 
new model because of the enhanced tendency of the giant 
planet to scatter material outward and to accrete material at 
late times during the migration. Accelerated protoplanetary 
accretion within the compacted shepherded portion of the 
disk is observed, with objects growing to several m e , but 
their orbits typically destabilize at late times due to capture 
in eccentricity-pumping mean motion resonances, result- 
ing in eventual collision with the giant planet. However, 
due to the sensitivity of our results to levels of dissipation 
close to the central star at late times, and the potential in- 
fluence of type I migration and eccentricity damping that 
we have neglected, it is premature to rule out the possibil- 
ity that hot-Neptunes or hot-Earths can form and survive 
interior to Jovian planets during their migration to become 
hot-Jupiters. 

4. Whilst we predicted terrestrial planet formation in hot- 
Jupiter systems from the results of our previous model, the 
external disks generated by our new model appear to be 
even more benign places for this to occur. These disks con- 
tain a greater mass of material, are less dispersed, and are 
composed of a higher fraction of small bodies capable of 
exerting dynamical friction and hence damping the excited 
orbits of scattered protoplanets. The formation of terrestrial 
planets of masses in the range 1 < m p < 3 m ffi occurred in 
or near the habitable zone during a simulation that we con- 
tinued after the giant stopped migrating. The radial mixing 
of volatile-rich material from beyond the snowline means 
that terrestrial planets forming in the habitable zones of 
hot-Jupiter systems are likely to be "water-worlds", host- 
ing deep, global oceans. 



The results presented in this paper make a clear prediction 
that terrestrial planets will eventually be discovered in the hab- 
itable zones of hot-Jupiter systems. Such systems may be de- 
tectable by forthcoming missions such as KEPLER, DARWIN 
and TPF. 
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